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Csl 1. INTRODUCTION 



The dominant mechanism for the production of a SM Higgs boson at the LHC is gluon-gluon fusion 
through a heavy (top) quark loop. For this reason this channel has attracted a large amount of theoretical 
attention [1]. Recently, the total cross section has been calculated to NNLO in the strong coupling 
constant a<j (i.e. at order 0{a'^)) [2-6] and also contributions from multiple soft gluon emission have 
been consistently included to NNLL accuracy [7]. In addition to the size of the total rate, a knowledge 
\ of the shape of the Higgs boson pT distribution is essential for any search and analysis strategies at the 

LHC. In particular, the pT distribution for the Higgs boson is expected to be harder than the one of its 



(N 



■ corresponding backgrounds. The Higgs boson px distribution has been computed with LL parton shower 

, Monte Carlos (HERWIG [8] and PYTHIA [9]), and through various resummed calculations. The latter 

^ \ techniques are the more powerful ones, but it is primarily the former that experimentalists at the LHC 

O ■ have to rely upon, because of their flexibility in allowing to test the effects of the various kinematic cuts 

which may optimize search strategies. 

In the kinematic region pj, <^ where most of the events are expected, large logarithmic 
Qh! corrections appear of the form a" In"* m|^/p|. that spoil the validity of the fixed order perturbative 

^ ' expansion. The pT distribution can be written as 

da da'"''- da^'^- 

+ ^^- (1) 



^ ■ dpj: dpj< dp. 

H ; 

The first term contains all logarithmically-enhanced contributions and requires their resummation to all 
orders. The second term is free from logarithmically-enhanced contributions and can be evaluated at 
fixed order in perturbation theory. The method to perform the all-order resummation is well known: to 
correctly take into account momentum conservation the resummation must be performed in the impact 
parameter (6) space [10, 11]. The large logarithmic contributions are exponentiated in the Sudakov form 
factor, which in the CSS [12] approach takes the form 
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Ac(a.(/i))ln (^j +Bc{as{p)) 



(2) 



where = 2e and c = q,g. The Ac and Be functions are free of large logarithmic corrections and 
can be computed as expansions in the strong coupling constant Ug : 
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Bc{as)= f:(-yB(-\ (4) 
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The functions Ac and Be control soft and flavour-conserving coUinear radiation at scales l/b^ fi^ rriH- 
Purely soft radiation at a very low scales /x^ 1/6 cancels out because the cross section is infrared safe 
and only purely collinear radiation up a scale jj, 1/6 remains, which is taken into account by the 
coefficients 

Cat{as,z) = Y,(^) Cj^\as,z). (5) 

n=l \ ^ / 

Beyond NLL accuracy, to preserve the process independence of the resummation formula, an additional 
(process dependent) coefficient H is needed [13], which accounts for hard virtual corrections and has an 
expansion 

F,(a,) = l + 5:(^] (6) 

n=l 

In the case of Higgs boson production through gg fusion, the relevant coefficients A^g\ A^g'^ and Bg^^ are 
known [14] and control the resummation up to NLL accuracy The NNLL coefficients C^f^ and Hg^^ 

(2) (3) 

are also known [13, 16]. The NNLL coefficient Bg has been computed in Refs. [17, 18], whereas Ag 
is not yet known exactly. In the following we assume that its value is the same that appears in threshold 
resummation [19]. 



2. PREDICTIONS FOR pr SPECTRA AND COMPARISONS 

In the 1999 Les Houches workshop, a comparison [20,21] of the HERWIG and PYTHIA (2 versions) 
predictions for the Higgs boson pT distribution with those of a pT resummation program (ResBos [22, 
23]) was carried out. This comparison was continued in the 2001 workshop and examined the impact 
of the B^'^^ coefficient [1]. In the meantime, a number of new theoretical predictions have become 
available, both from resummation and from the interface of NLO calculations with parton shower Monte 
Carlos. For these proceedings, we have carried out a comparison of most of the available predictions for 
the Higgs boson px distribution at the LHC. We have used a Higgs boson mass of 125 GeV and either 
the MRST2001 or the CTEQ5M pdf 's. The difference between the two pdf 's for the production of a 125 
GeV mass Higgs boson is of the order of a few percent. Before comparing the different predictions, we 
comment on the various approaches in turn. 

Parton shower MC programs such as HERWIG, which implements angular ordering exactly, im- 
plicitly include the A^^\ A^"^^ and B^^^ coefficients and thus correctly sum the LL and part of the N'^LL 
contributions. However, in the most straightforward implementations, MC cannot correctly treat hard 
radiation. By contrast, the PYTHIA MC, which does not provide an exact implementation of angular 
ordering, has a hard matrix element correction ^. Recently, an approach to match NLO calculations to 
parton showers generators, MC@NLO [25, 26], has been proposed, and applied, amongst the other, to 
Higgs production. This method joins the virtues of NLO parton level generators (correct treatment of 
hard radiation, exact NLO normalization) to the ones of MC. It thus can be compared to a resummed 
calculation at NLL-i-NLO accuracy. 

As far as resummed calculations are concerned, we first consider two implementation of the CSS 

(12 3) fl2) (1) 

approach. The ResBos code includes the Ac , Be and C^^ coefficients in the low-px region and 
matches this to the NLO distribution at high px- NNLO effects at high px are approximately taken into 
account by scaling the second term in Eq. Q with a K-factor. The matching is performed through a 
switching procedure whose uncertainty will be considered in the following. The calculation of Berger 
and Qiu [27] also performs a pT resummation in b space and is accurate to NLL. The coefficient B^'^^ 
is included but the matching is still to NLO. Note that in both these approaches the integral of the 

'There are two different classification schemes of the LL, NLL, NNLL, etc terms and their corresponding B contents. Here 
we use the most popular scheme. Another is discussed in Ref. [15]. 

^Very recently hard matrix element corrections for Higgs productions have been implemented in HERWIG as well [24]. 



spectrum is affected by higher-order contributions included in a non-systematic manner whose effect is 
not negUgible for Higgs production. 

The prediction by Bozzi, Catani, de Florian and Grazzini [28] (labeled Grazzini et al. in the 
following) is based on an implementation of the 6-space formalism described in [13,28]. The calculation 
has the highest nominal accuracy since it matches NNLL resummation at small pT to the NNLO result 

(2) (2) 

at high pt [29]. This approach includes the coefficients C^^ and Hg in approximated form. The 
main differences with respect to the standard CSS approach are the following. A unitarity constraint is 
imposed, such that the total cross section at the nominal (NNLO) accuracy is exactly recovered upon 
integration. A study of uncertainties from missing higher order contributions can be performed as it is 
normally done in fixed order calculations, that is, by varying renormalization and factorization scales 
around the central value, that is chosen to be mn- 

Finally, we discuss the px distribution of Ref. [30] (Kulesza et al.). This is obtained using a joint 
resummation formalism, by which both threshold and low-pr logarithmic contributions are resummed 
to all orders. This approach has been formally developed to NLL accuracy, but the NNLL coefficients 
A(3),i?(2),cW and i^W can also be incorporated. The matching is still performed to NLO. Even 
though a low mass Higgs boson at the LHC is produced with relatively low x partons, threshold effects 
can still be significant due to the large color charge in the gg initial state as well as steep x dependence 
of the gluon distribution functions at low x. This leads to an increased sensitivity to Sudakov logarithms 
associated with partonic threshold for gluon-induced processes, as shown in Ref. [7]. 

It is known that the low-px region is sensitive to non-perturbative effects. These aie expected 
to be less important in the gluon channel due to the larger colour charge of the gg initial state [21]. 
Different treatments of non-perturbative effects are included in the ResBos, Berger et al. and Kulesza et 
al calculation, whereas Grazzini et al. prediction is purely perturbative. 

The absolute predictions for the cross sections are shown in Figure ^ All curves are obtained in 
the m-top oo limit. HERWIG and PYTHIA cross sections are significantly smaller than the other 
predictions, their normalization being fixed to LO. In the high-pr region, the HERWIG prediction drops 
quickly due to the lack of hard matrix element coiTcctions. PYTHIA, in contrast, features the hard matrix 
element corrections. We also note that PYTHIA prediction is significantly softer than all the other curves, 
and thus its overall shape is fairly different from all the other predictions. 

The MC@NLO cross section, about 32.4 pb, is roughly twice that of the HERWIG and PYTHIA 
predictions, being fixed to the NLO total cross section. 

Two predictions (step, smooth) are shown for ResBos which differ in the manner in which the 
matching at high px is performed. Their difference can be considered as an estimate of the ambiguity 
in the switching procedure. The two curves correspond to the same total cross section of about 36.2 
pb, which is about 8 % higher than the NLO cross section. This is the effect of the higher-order terms 
that enter the prediction for the total rate in the context of the CSS approach. A slightly softer curve is 
obtained by Berger and Qiu. The predicted cross section (37 pb) is close to that of ResBos. 

The Grazzini et al. prediction has an integral of about 39.4 pb, which corresponds to the total 
cross section at NNLO. Contrary to what is done in Ref. [28], here the curve is obtained with MRST2002 
NNLO partons and three-loop a^. The difference with the result obtained with MRST2001 NNLO PDFs 
is completely negligible. 

Concerning the Kulesza et al curve, the subleading terms associated with low x emission (i.e. in 
the limit opposite to partonic threshold) and of which only a subset is included in the joint resummation 
formalism, play an important role numerically. As a result, the total cross section turns out to be 35 pb, 
about 10% lower than the pure threshold result, which is 39.4 pb [30]. 

We now want to examine in more detail the relative shapes of the predictions plotted in Figure. Q 
In Figure. IJlall the predictions are normalized to the Grazzini et al. cross section of 39.4 pb. 

In the region of small and moderate pT (say, smaller than 100 GeV) all of the predictions are 



basically consistent with each other, with the notable exception of PYTHIA, which predicts a much 
softer spectrum. The curve of Kulesza et al. is also softer than the others. 

For larger pT, HERWIG gives unreUable predictions, since the transverse momentum is generated 

solely by means of the parton shower, and therefore it lacks hard matrix element effects. The Grazzini et 
al. and ResBos curves are harder than MC@NLO for large pT- There are two reasons for this. Grazzini 
et al. implement the NNLO matrix elements exactly, corresponding to the emission of two real partons 
accompanying the Higgs in the final state [29]; ResBos mimics these contributions, by multiplying the 
NLO matrix elements by the K factor. MC@NLO, on the other hand, contains only NLO matrix elements 
(one real parton in the final state). Secondly, Grazzini et al. and ResBos choose the renormalization and 
factorization scales equal to mu, whereas in MC@NLO these scales are set equal to the transverse mass 
of the Higgs, ^ rnj^ + pj<. The difference is small at the level of total rates, but it is not negligible in the 
tail of the pt distribution. 

3. CONCLUSIONS 

Up to now, the ATLAS and CMS experiments have relied primarily on the predictions from HERWIG 
and PYTHIA in designing both their experiments as well as defining their search and analysis strate- 
gies. In the last few years, a number of tools for and predictions of the Higgs boson cross section at the 
LHC have become available, with the inclusion of beyond-the-leading-order effects at different level of 
accuracy. In the case of total rates, NNLO results have recently become available; their consistent inclu- 
sion in experimental analysis will allow to further decrease the estimated lower bound on the integrated 
luminosity to be collected for discovery. 

In this contribution, we primarily focused on the predictions for the pT spectrum, comparing the 
results of Monte Carlos with those obtained with analytically-resummed calculations. In contrast to 
the situation in 1999, all of the predictions, with the exception of PYTHIA, result in the same general 
features, most notably in the position of the peak. However, differences do arise, because of different 
treatments of the higher orders. It is an interesting question beyond the scope of this review that of 
whether these differences are resolvable at the experimental level, which may lead to modify the strategy 
for searches. In order to answer this, studies including realistic experimental cuts must be performed 
with the newly available tools. 

References 

[I] D. Cavalli et al. 2002. 

[2] S. Catani, D. de Florian, and M. Grazzini. JHEP, 05:025, 2001. 

[3] R.V. Harlander and W.B. Kilgore. Phys. Rev., D64:0I30I5, 2001. 

[4] R.V. Harlander and W.B. Kilgore. Phys. Rev. Lett, 88:201801, 2002. 

[5] C. Anastasiou and K. Melnikov. Nucl. Phys., B646:220-256, 2002. 

[6] V. Ravindran, J. Smith, and W. L. van Neerven. Nucl. Phys., B665:325-366, 2003. 

[7] S. Catani, D. de Florian, M. Grazzini, and R Nason. JHEP, 07:028, 2003. 

[8] G. Corcella et al. 2002. 

[9] T. Sjosti-and, L. Lonnblad, S. Mrenna, and R Skands. 2003. 

[10] G. Parisi and R. Petronzio. Nucl. Phys., B154:427, 1979. 



[11] 



YL. Dokshitzer, D. Diakonov, and S. I. Troian. Phys. Rept, 58:269-395, 1980. 



[12] J.C. Collins, D.E. Soper, and G. Sterman. Nucl. Phys., B250:199, 1985. 
[13] S. Catani, D. de Florian, and M. Grazzini. Nucl. Phys., B596:299-312, 2001. 
[14] S. Catani, E. D'EmiUo, and L. Trentadue. Phys. Lett., B21 1:335-342, 1988. 
[15] W. Giele et al. 2002. 

[16] R. P. Kauffman. Phys. Rev., D45:1512-1517, 1992. 

[17] D. de Florian and M. Grazzini. Phys. Rev. Lett, 85:4678^681, 2000. 

[18] D. de Florian and M. Grazzini. Nucl. Phys., B616:247-285, 2001. 

[19] A. Vogt. Phys. Lett., B497:228-234, 2001. 

[20] S. Catani et al. 2000. 

[21] C. Balazs, J. Huston, and I. Puljak. Phys. Rev., D63:014021, 2001. 

[22] C. Balazs and C. P. Yuan. Phys. Rev., D56:5558-5583, 1997. 

[23] C. Balazs and C. P Yuan. Phys. Lett., B478: 192-198, 2000. 

[24] G. Corcella and S. Moretti. 2004. 

[25] S. Frixione and Bryan R. Webber. JHEP, 06:029, 2002. 

[26] S. Frixione, R Nason, and B.R. Webber. JHEP, 08:007, 2003. 

[27] E.L. Berger and J. Qiu. Phys. Rev., D67:034026, 2003. 

[28] G. Bozzi, S. Catani, D. de Florian, and M. Grazzini. Phys. Lett., B564:65-72, 2003. 
[29] D. de Florian, M. Grazzini, and Z. Kunszt. Phys. Rev. Lett, 82:5209-5212, 1999. 
[30] A. Kulesza, G. Sterman, and W. Vogelsang. Phys. Rev., D69:014012, 2004. 





20 40 60 80 100 120 140 160 180 jq^^OO 

Fig. 1: The absolute predictions for tlie production of a 125 GeV mass Higgs boson at the LHC. 



